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ABSTRACT 

Thermodynamics of helix-coil transitions in amino-acid homo-oligomers are studied by 
the recently proposed multicanonical algorithms. Homo-oligomers of length 10 are consid- 
ered for three characteristic amino acids, alanine (helix former), valine (helix indifferent), 
and glycine (helix breaker). For alanine other lengths (15 and 20) are also considered 
in order to examine the length dependence. From one multicanonical production run 
with completely random initial conformations, we have obtained the lowest-energy con- 
formations and various thermodynamic quantities (average helicity, Zimm-Bragg s and 
a parameters, free energy differences between helix and coil states, etc.) as functions of 
temperature. The results confirm the fact that alanine is helix-forming, valine is helix- 
indifferent, and glycine is helix-breaking. 
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INTRODUCTION 

Recent experimental measurements (for a review, see Ref. of the a-helix propensities 
of amino acids in short peptide systems have raised a renewed interest in the theoretical 
studies of a-helix formation, and a number of simulation results have been reported. M- 
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The major difficulty in conventional protein simulations such as molecular dynamics 
lies in the fact that simulations at temperatures of experimental interest tend to get 
trapped in one of a huge number of local minima of potential energy surface. Hence, 
the simulations strongly depend on the initial conditions. This is why most simulations 
start from a folded conformation that is suggested by X-ray or NMR experiments and the 
unfolding of the conformation is studied. However, this is a serious limitation if one is 
interested in the prediction of protein structures from the first principles without the use of 
information on experimentally implied structure. Hence, novel algorithms that accelerate 
thermalization are in urgent demand. Recently, the authors proposed the application of 
the multicanonical algorithms [|13|, [Uj to the protein folding problem. |15|] The performance 
of the algorithm was compared with that of Monte Carlo simulated annealing, ||16|| another 
effective method for overcoming the above-mentioned multiple-minima problem, and it 
was claimed that the former is superior to the latter. [[Ij], [H| The same algorithm was 



referred to as entropy sampling by another group, [T^] but the proof of the equivalence 
of the two methods was given to clarify the matter. [EDI Apart from the protein folding 
problem the multicanonical approach was also successfully applied to the similar problem 



of spinglasses. |21||- [[23|| The advantage of this new algorithm lies in the fact that it not 
only alleviates the multiple-minima problem but also allows the calculation of various 
thermodynamic quantities as functions of temperature from one simulation run. The 
purpose of the present work is to further test the effectiveness of the new algorithm in 
the study of thermodynamics of the protein folding problem. 

In this article, we study thermodynamics of helix-coil transitions in amino-acid homo- 



oligomers by multicanonical algorithms. Preliminary results were reported elsewhere. [24 



Homo-oligomers of length 10 are considered for three characteristic amino acids, alanine 
(helix former), valine (helix indifferent), and glycine (helix breaker). We first investigate 
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the lowest-energy conformations obtained by our simulations. We then calculate various 
thermodynamic quantities (such as the average % helix, Zimm-Bragg s and a parameters, 
free energy differences of helix-coil transitions, etc.) over a wide range of temperatures 
(from 100 K to 1000 K). To our knowledge this is the first time that such a wide range 
of temperatures could be covered by a single simulation run to calculate various ther- 
modynamic quantities. The Zimm-Bragg parameters are compared with that of recent 
experiments. 
METHODS 

Multicanonical Ensemble 

Although the algorithms are explained in detail elsewhere (see, for instance Refs. [ |I4| , |T5| , 
18]]), we briefly summarize the idea and implementations of the method for completeness. 

Most Monte Carlo (MC) simulations are done in the canonical ensemble which is 
characterized by the Boltzmann weight factor wb(E) = exp (^—(3E^j . Here (3 = 1/RT 
is the inverse temperature. States of energy E are then distributed according to the 
probability 

P B (T,E) oc n(E)w B (E) , (1) 

where n(E) is the density of states. Since n(E) is a rapidly increasing function and the 
Boltzmann factor decreases exponentially, Pb{T,E) generally has a bell-like shape with 
its value varying many orders of magnitude as a function of E. On the other hand, in 
the multicanonical approach Monte Carlo simulations are performed in an artificial multi- 
canonical ensemble, [13|] which is definedby the condition that the probability distribution 



of the energy shall be constant: 

P mu (E) oc n(E)w mu (E) = const. (2) 

All energies have equal weight and a one-dimensional random walk in energy space is 
realized, which insures that the system can overcome any energy barrier. Note that from 
Eq. (0) we have 

w mu {E) oc n-\E) . (3) 

Unlike in the case for the canonical ensemble, the multicanonical weight factor w mu (E) 
is not a priori known, and one needs its estimator for a numerical simulation. Hence, 



the multicanonical ansatz consists of three steps: In the first step the estimator of the 
multicanonical weight factor w mu (E) is calculated. Then one performs with this weight 
factor a multicanonical simulation with high statistics. The standard Markov process (for 
instance, in a Metropolis update scheme ||25|| ) is well-suited for generating configurations 
which are in equilibrium with respect to the multicanonical distribution. Monitoring the 
energy in this simulation one would see that a random walk between high energy states 
and ground-state configurations is realized. In this way information is collected over the 
whole energy range. Finally, from this simulation one can not only locate the energy 
global minimum but also obtain the canonical distribution at any inverse temperature j3 
for a wide range of temperatures by the re- weighting techniques: [f?Gf 

P B (T,E)<xP mu (E)w^ u e^ E . (4) 

This allows one to calculate any thermodynamic quantity at temperature T. For instance, 
the expectation value of a physical quantity O at temperature T is given by 



dE 0(E)P B (T,E) 

<0> T = J — r . (5) 

J dE P B (T,E) 

In the ideal case a re-weighting is possible to all temperatures. However, in reality it may 
not be possible or useful to ensure the condition of Eq. (H) for all energies, but only for 
an interval E min < E < E max . In this case, the range of temperatures 

T^min — T — T max (6) 

for which the re-weighting yields correct expectation values has to be determined from 



the condition HT4 



E m in < < E > T < E max . (7) 

The crucial point is the first step: calculating the estimator for the multicanonical 
weight factor w mu (E). In Ref. |1§ we proposed for this purpose the following iterative 



procedure: 

1. Perform a canonical Monte Carlo simulation at a sufficiently high temperature To. 
In our case we chose T = 1000 K. The weight factor for this simulation is given 
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by w B (E) = e-P° E with (3 = 1/BT . Initialize the array S(E) to zero, where E is 
discretized with bin width SE (= 1 kcal/mol in the present work). 

2. Sample the energy distribution obtained in the previous simulation as a histogram 
H(E) with the same bin width as for S(E). In the first iteration (step 1 above) 
determine E max as the value near the mode where the histogram has its maximum 
{E-max is fixed throughout the iterations). Let E min be the lowest energy obtained 
throughout the preceding iterations. For all H(E) with entries greater than a certain 
minimum value (say, 20) and E min < E < E max , update the array S(E) by 

S{E) = S{E) + In H (E) . (8) 

3. Calculate the following multicanonical parameters a(E) and (3{E) from the array 
S(E): 

A) > E > E max 

P(E) = \ ft + S{E 'l - S J E) , E mm < E < E> < E max (9) 



E' — E 

(3(E min ) , E < E,, 



mm 



and 



/ 7ti\ J , E ^ E max (~\n\ 

a ^> - \ a(E') + {P{E') -P(E))E' , E<E n ' ' 



■'max 



where E and E' are adjacent bins in the array S(E). 

4. Start a new simulation with the multicanonical weight factor defined by 

w mu {E) = e-^ E ~^ . (11) 

5. Iterate the last three steps until the obtained distribution H(E) becomes reasonably 
flat in the chosen energy range. 

While this method for determining the multicanonical weight factor w mu (E) is quite 
general, it has the disadvantage that it requires a certain number of iterations which is 
not a priori known. For the calculations in Ref. and Ref. [I8|, about 40 % and 4 
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of the total CPU time was respectively spent for this task. We remark that the above 
method of calculating multicanonical weights is by no means unique. Especially it is not 
necessary to choose the parametrization of Eq. (|TT| ) for the multicanonical weight factor. 



However, with this parametrization and its introduction of "effective" temperatures (3(E) 
the connection to the canonical ensemble becomes apparent. Especially, if the parameters 
(3(E) of Eq. (0) are a monotone function of the energy E, then they indicate the range 
of temperatures, for which a valid re- weighting is possible (see Eq. (|7|)). Note also that 
S(E) for (3o = is an estimator of the microcanonical entropy. 
Multicanonical Annealing 

If one is just interested in the ground-state structure, it may be worthwhile to use instead 
a variant of the multicanonical method, multicanonical annealing. |27|, |17], |18| Here, an 
upper bound in energy is introduced at the other end of the annealing direction, rejecting 
all attempts beyond this bound. Annealing is achieved by moving the bound in the 
annealing direction while keeping the sampling interval AE fixed. Within this energy 
interval the system can move out of local minima as long as their barrier heights do not 
exceed the upper limit of the energy range. Because of the finite interval size AE, the 
MC procedure will no longer be ergodic and it is not possible to find the equilibrium 
properties of the system. Hence, unlike in the case for regular multicanonical algorithm, 
the canonical distribution cannot be reconstructed. For the purpose of annealing this does 
not matter as long as one chooses the sampling interval large enough to allow important 
fluctuations throughout the annealing process. However, since ergodicity is not fulfilled, 
one has to repeat the annealing process several times with different initial conformations 
to make sure that one has found a good approximation to the global minimum. The 
optimal sampling interval AE is not known a priori for multicanonical annealing and has 
to be chosen on a trial-and-error basis. 



The multicanonical annealing algorithm is discussed in detail in Refs. JT7, [18] where 
first tests of the method for the protein folding problem were performed. Results better 
than those obtained by simulated annealing were reported. The following implementation 
of the algorithm was proposed: 

1. Perform a short canonical Monte Carlo simulation at a sufficiently high temperature 
To. Again we chose To = 1000 K in the present work. Initialize an array S(E) to 
zero, where E is discretized with bin width SE (= 1 kcal/mol in the present work). 
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2. Sample the energy distribution obtained in the previous simulation as a histogram 
H(E). Let E min be the lowest energy obtained throughout the preceding iterations. 
For all H(E) with entries greater than a certain minimum value (say, 20), update 
the array S(E) by 

S(E) = S{E) + lnH{E) . (12) 

3. Update the upper bound E wa u of the sampling interval by 

E wa ii = max(Ei ast , E min + AE) , (13) 

where AE is the size of the sampling energy range and Ei ast the energy of the last 
conformation. 

4. Calculate the following parameter /3 min by 

R _ S{E waU ) - S(E min ) 

Pmin ~ rp _ rp " I 1 V 

£-<watt -E-'min 

5. Define the new weight factor by 

{0 , E > E wa u 

e~ s{E) _ , E mm <E<E wall (15) 

c j J- 1 ^ J -'min 

and perform a new simulation with this weight factor, starting from the last con- 
formation of the preceding simulation. 

6. Iterate the last four steps till no newer E m i n is found for a certain number of con- 
secutive iterations. 

Peptide Preparation 

We considered amino-acid homo-oligomers of Ala, Val, and Gly. By experiments Ala is 
known to be a strong helix former, while Val and Gly are known to be helix indifferent 
and helix breaker, respectively. The number of residues, N, for each homo-oligomer was 
taken to be 10. For (Ala)^, however, the cases for N = 15 and 20 were also considered 
in order to examine the N dependence. Since the charges at peptide termini are known 
to reduce helix content, p8| p9fl we removed them by taking a neutral NH2- group at the 
N-terminus and a neutral -COOH group at the C-terminus. 



Potential Energy Function 

The potential energy function E tot that we used is given by the sum of the electrostatic 
term Ec, 12-6 Lennard- Jones term Elj, and hydrogen-bond term Ehb for all pairs of 
atoms in the peptide together with the torsion term E tor for all torsion angles: 

E to t = E c + E LJ + E HB + E tor , (16) 
Ec = (17) 

= E(p£-fO. (18) 

^or = E^( 1±C0S ( n ^/))- (20) 

I 

Here, is the distance between the atoms i and j, and xz is the torsion angle for 
the chemical bond /. The parameters (q i7 Aij, B^, G^, D^, Ui, and n;) for the energy 
function were adopted from ECEPP/2.|3(|-[[3^] Since one can avoid the complications of 
electrostatic and hydrogen-bond interactions of side chains with the solvent for nonpolar 
amino acids, explicit solvent molecules were neglected and the dielectric constant e was 
set equal to 2. This is a crude approximation, and we have to keep this limitation in mind 
when we compare our simulation results with experiments. The computer code KONF90 
was modified to accommodate the multicanonical algorithms. (There are slight 



differences in conventions between KONF90 and the original version of ECEPP/2; for 
example, of KONF90 is equal to <p\ — 180° of ECEPP/2, and energies are also different 
by small irrelevant constant terms.) The peptide-bond dihedral angles uj were fixed at 
the value 180° for simplicity, which leaves fa, ipi, and Xi {i — 1> ' ' " > N) as independent 
degrees of freedom. Since Ala has one x i n the side chain, Val three x' s > an d Gly none, 
the numbers of independent degrees of freedom are 30, 50, 20, 45, 60 for (Ala)io, (Val) 10 , 
(Gly) io, (Ala) 15, (Ala) 20 , respectively. 
Computational Details 

One MC sweep updates every dihedral angle (in both the backbone and the side chains) 
of the homo-oligomers once. 

For regular multicanonical simulations, the multicanonical weight factors were deter- 
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mined by the iterative procedure described above. We needed between 40,000 sweeps 
(for (Ala)io) an d 100,000 sweeps (for (Val) 10 ) for their calculation. All thermodynamic 
quantities were then calculated from one production run of 200,000 MC sweeps following 
additional 10,000 sweeps for equilibration. 

For multicanonical annealing, we performed 10 simulation runs from different random 
initial conformations with 20,000 MC sweeps (so that the total number of MC sweeps is 
equal to 200,000 for each homo-oligomer) . We divided the 20,000 MC sweeps of each run 
in 10 annealing iterations of 2,000 MC sweeps. The value of the sampling interval AE 
was chosen to be AE = 15 kcal/mol. 

In all cases, each simulation started from a completely random initial conformation 
("Hot Start"). We remark that for the regular multicanonical algorithm, simulations with 
initial conformations of an ideal helix ("Cold Start") were also performed, and we found 
that the results are in agreement with those from random initial conformations. This 
suggests that thermal equilibrium has been attained in our simulations. 
RESULTS AND DISCUSSION 
Lowest-Energy Conformations 

We first investigate the lowest-energy conformations obtained from our simulations. Here, 
we are able to cross-check our results of regular multicanonical simulation by those of 



multicanonical annealing. The criterion we adopt for a-helix formation is as follows: [11 



We consider that a residue is in the right-handed a-helix configuration when the dihedral 
angles ((f), ip) fall in the range (—70 ± 20°, —37 ± 20°). The length in of a helical segment 
is then defined by the number of successive residues which are in the right-handed a-helix 
configuration. The number tir of helical residues in a conformation is defined by the 
sum of £r over all right-handed helical segments in the conformation. Note that we have 
n R > with the equality holding when there exists just one (or no) helical segment in 
the conformation. Since Gly has no side chain, it can produce both right-handed and 
left-handed helices. Analogous quantities £l and ul for left-handed a-helix are defined 
with obvious reversal of signs for the dihedral angles. 

The conformations obtained during the simulation are classified into two states, helix 
and coil. Here, a conformation is considered to be in the helix state, if it has a segment 
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with helix length t > 3. A conformation is considered to be in the coil state, if it is not 
in the helix state. We remark that this definition of the helical conformation is in some 
sense arbitrary; for instance, we could define the helical conformation as that with I > 2 
instead. But one has to draw a line somewhere, and we chose this definition, since i = 3 
corresponds to roughly one turn of a-helix. In a similar way, the above definition for the 
helical residue is also by no means unique. We have chosen the central values of —70° 
and —37° here, because they were the average values of (j) and ip in the helical residues 



we obtained with KONF90 in previous simulations. We have to keep in mind these 
ambiguities in the definitions of helical residue and conformation, but we remark that we 
found no qualitative changes when we checked our results by using different definitions in 
our analysis. 

In Table 1 we list the energy E (kcal/mol) and the maximum helix lengths £r and £l of 
the lowest-energy conformations obtained during each of 10 multicanonical annealing runs 
of 20,000 MC sweeps for (Ala)i , (Val)i , and (Gly)i . Each run was performed from a 
completely random initial conformation. Hence, the results of the 10 runs are independent 
of each other. These quantities are listed separately for conformations in both helix and 
coil states. Our aim was to study if there is a unique ground state or a multitude of low- 
energy states. Despite the experimentally observed differences in a-helix propensities, 
one could assume that an ideal helix is the unique ground-state structure for all homo- 
oligomers, since this structure is energetically favored by the van der Waals term. The 
observed propensity differences would then be caused by the entropic effects. To check this 
conjecture we also calculated the minimized energies of an idealized helical conformations 
for the three homo-oligomers and compared them with our numerical results. To obtain 
the idealized helical structure we first set all backbone angles to <fi = —70° and ip = —37° 
(i.e., £ = 10) and the side chain angles x to 60° for (Ala)io and 180° for (Val)io- The 
energies of these initial conformations were then minimized by the Newton method. The 
resulting energies are also listed in Table 1. In Table 2 we give the dihedral angles of the 
lowest-energy conformations of Table 1 for completeness. 

The first thing one can tell from Table 1 is that the lowest-energy structures of (Ala) io 
are all right-handed a-helix with almost 100 % helicity (£r = 8). The energy of these 
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conformations is comparable to that of the minimized ideal helical conformation. Note 
that the energy of the minimized ideal-helix structure is slightly higher than the lowest 
energy of our multicanonical annealing simulation. We conclude that the structure with 
the global- minimum energy for (Ala)io is an ideal helix and that there is a continuum 
of excited states whose energy and structure are similar to those of the ground state. 
We remark that with the present multicanonical annealing simulation, the probability 
of finding the ground-state structure (£r > 8) is 60 %; 6 out of 10 runs found it. The 
probability would be higher if we increased the statistics of each run. Finally, we found 
that the side chain structure of (Ala)io is also unique for the ground-state structure; 
namely, the values of x are one of 60°, —60°, and 180°, which are all equivalent angles 
because of the three-fold rotational symmetry of the alanine side chain (see Table 2). 

Things look different for the other two homo-oligomers. Since Gly has no side chain, 
it can produce both right-handed and left-handed helices. The longest left-handed helical 
length £ L for (Gly) io is as large as 6 in Table 1. On the other hand, the longest right- 
handed helical length in for (Gly)io is only 3 in Table 1. Here, it appears that with 
our energy function, left-handed helix is favored over right-handed one. However, as 
is discussed below, the regular multicanonical production run did produce longer right- 
handed helix {Ir = 5) with similar energy (see Table 3 below). Hence, this tendency 
is weak and caused by the lack of statistics in our multicanonical annealing simulations. 
The energy of the lowest-energy helical conformation is higher than or comparable with 
that of coil structures. The lowest-energy structure is a coil (run 8). Furthermore, the 
energy of the minimized ideal-helix structure is also higher than that of the lowest-energy 



coil structure. Similar observations were made in Ref. |35| and it was conjectured that 
in a short a-helix the electrostatic interaction is unfavorable due to parallel arrangement 



of peptide bond dipoles. |35| However, it was argued there that for larger chains, the 
attractive van der Waals term will win and a helix conformation would be the ground 
state. From our multicanonical annealing results we can only conclude that for (Gly)io 
a-helix is not the ground-state conformation and that there exist a multitude of different 
low-energy (coil) states. Without side chains, (Gly)io is so flexible that the ground state 
cannot be strongly energetically favored compared to other structures. In this sense 
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(Gly)io is like a spinglass with many local minima of energy near the global minimum. 

For (Val) 10, some of the low-energy structures are helical and others not. Our lowest- 
energy conformation (run 9) is the one with the highest value of £r (= 6). As is shown 
in Table 1, the lowest-energy conformations obtained by the multicanonical annealing 
runs have a much higher energy (at least 6 kcal/mol) than the minimized ideal-helix 
structure. Hence, we have to admit that the multicanonical annealing simulation did not 
find the energy global minimum for (Val)io, although it got close to it. We conjecture 
that for (Val) io an ideal helix structure is the ground state and that because of the steric 
hindrance of the side chains (Val has a large side chain), it is separated by high energy 
barriers from a multitude of states with smaller helical lengths and slightly higher energy. 
Hence, it is very difficult to reach the ideal helix structure from a completely random 
initial conformation. 

We now compare the above results by multicanonical annealing with those by regular 
multicanonical algorithm production runs. For each homo-oligomer, (Ala)i , (Ala) 15, 
(Ala) 20, (Val)io, and (Gly)io, one production run with 200,000 MC sweeps (after 10,000 
sweeps for equilibration) was made. Hence, the statistics is 10 times more than one 
multicanonical annealing run (but 10 runs were made for the latter). In Table 3 we list 
the energy E (kcal/mol), helix lengths £r and £l, and dihedral angles of the lowest-energy 
conformations obtained by each multicanonical production run. These quantities are again 
listed separately for the lowest-energy conformations in both helix and coil states. As is 
clear from the table, the uniqueness of the global-minimum state (ideal helix structure) 
for (Ala) 10 is apparent in accord with the implication of the multicanonical annealing 
runs; we have roughly the same energy (es —10 kcal/mol), helix length (£r = 8), and 
dihedral angles in Tables 2 and 3. We find that even the structure of the termini is 
unique; the dihedral angles ipi, 0io, and ^10 are essentially the same in Tables 2 and 
3. We remark that the lowest-energy structures for (Ala) 15 and (Ala) 2 o obtained by the 
multicanonical production runs were also ideal helix with helix lengths £r — 13 and 
18, respectively. Furthermore, the values of x were one of 60°, —60°, and 180°. Non- 
uniqueness of the lowest-energy structures for (Val) 10 and (Gly)io are also clear in the 
table; we have different helix length and the dihedral angles but similar energies (~ 2, —8 
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kcal/mol for (Val)io, (Gly)io, respectively) in Tables 2 and 3. Again, the global-minimum 
conformation for (Val) 10 (ideal helix) was not reached. However, we remark that the 
lowest energies for helix states reached by multicanonical production runs for (Ala)io and 
(Gly)io are less than the minimized ones in Table 1. In Figure 1 we show stereoscopic views 
of the lowest-energy structures from Tables 2 and 3 for completeness. The lowest-energy 
left-handed a-helix structure for (Gly)io (from Table 3) is also shown in the figure. 

As one can read off from Table 3, the differences AE = Ec — Eh of the lowest energies 
between coil state (C) and helix state (H) are 12.4, 4.8, —0.9 kcal/mol for (Ala)io, (Val) 10 , 
(Gly)io, respectively. (Note that these values are lower bounds and the more realistic 
values may be even larger, since their coil states in Table 3 are almost helical. Namely, if 
we define a helical residue by the condition (<p,ip) = (—70 ±23°, —37 ±23°) instead, these 
structures are considered to be helical.) This and our results from the multicanonical 
annealing runs imply that the helix state is energetically favored strongly for (Ala)io and 
slightly favored for (Val) 10 , while it is not favored for (Gly)io.[|ll| These observed energy 
differences are a natural explanation for the experimentally observed fact that Ala is 
a helix former and Gly is a helix breaker, while Val comes in between the two. More 
quantitative arguments about this point will be given in the following subsections. 
Energy Distribution and Specific Heat 

While multicanonical annealing gives information only on the lowest-energy state, regular 
multicanonical algorithm also allows the calculation of thermodynamic quantities at vari- 
ous temperatures. The range of temperatures, for which a valid re-weighting is possible, is 
given in Eq. (^). Since we are fairly sure that we determined the ground state for (Ala)io, 
we can in principle trust our results here down to T = K. However, we do not know 
for sure that we found the ground state for (Gly)io and we know that we missed it for 
(Val) io- Hence, we had to determine in these cases the range of temperatures that allow 
reliable re- weighting by Eq. (|7]). As an upper limit to which we can trust our results for 
the three homo-oligomers, we found T m , n = 200 K. Therefore we restricted our analysis 
to temperatures between 200 K and 1000 K in most of the cases below. 

In Figure 2 we show 3-dimensional plots of the probability distribution of energy as 
a function of temperature for (Ala)i , (Val)i , and (Gly)i . These results were obtained 
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from the same multicanonical production runs as above with the reweighting techniques of 
Eq. (4). At each fixed temperature T, the probability distribution P B (T,E) corresponds 
to a canonical distribution at this temperature. Thus, it is given by Eq. (1) and has a 
bell-shape. The width of the bell-shape is large for high temperatures reflecting the large 
energy fluctuations. As the temperature decreases, this width decreases but the height 
of the peak increases (since the probability is normalized to 1), and it should behave as 
a 5-function 5(E — Eq) in the limit T — > K, where Eq is the global-minimum energy. 
These properties of canonical ensemble are indeed clearly seen in Figure 2. 

In order to demonstrate the reliability and superiority of our method, we also per- 
formed for (Ala) io and (Ala) 2 o canonical simulations of 200,000 sweeps at temperature 
T = 270 K and compared the obtained distribution of energies with the one from multi- 
canonical simulation by reweighting. Again 10,000 sweeps were performed for equilibra- 
tion of our system. Figures 3 show our results. Even for N = 10 we observe remarkable 
differences in the shape of the distribution, and they are totally different for N = 20. In 
the case of N — 10 we observe for the canonical simulation a pronounced tail for energies 
larger than 1 kcal/mol while both distribution look similar for lower energies. When we 
looked deeper into the time series for the canonical run, we found that our system need 
at least 80,000 sweeps more than the estimated 10,000 sweeps for equilibration. When we 
discard these sweeps, the distribution agreed with that obtained by the multicanonical 
production run. In the case for N = 20 the system did not equilibrate at all in the canon- 
ical simulation. On the other hand, the multicanonical simulation was fully equilibrated: 
The distribution did not depend on the initial conformations. Moreover, our method gives 
information on a whole temperature range, not only a single one. This demonstrates the 
superiority of our method. 

In Figure 4 we show the "specific heat" as a function of temperature for the three 
homo-oligomers, (Ala)i , (Val)i , and (Gly)i . The specific heat here is defined by the 
following equation: 

C0)=P 2 <Elt> ~ N <EtOt> \ (21) 

where j3 is the inverse temperature 1/RT and N (= 10) is the number of residues in the 
oligomers. It was calculated from the reweighting techniques of Eqs. (4) and (5). The 
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data in the figure all represent a peak at a certain temperature, indicating that there is a 
crossover between a coil phase and a helix phase. The temperatures at the peak, transition 
temperatures, are T c « 430, 330, 360 K for (Ala)i , (Val)io, (Gly)io, respectively. The 
peak structure is most conspicuous for (Ala)io, an d this corresponds to a transition from 
a random-coil phase to an ideal-helix phase. The transition temperature T c for (Ala)io is 
rather high. This indicates that (Ala)io is substantially helical (see Figure 7 below) even 
at temperatures near 400 K. The peaks in specific heat are not as sharp for (Val)io and 
(Gly)io, and this is a typical characteristic of a spinglass. 

To study if the mentioned crossover is an indication for a "phase transition" we show in 
Figure 5 the specific heat as a function of temperature for (Ala) 10 , (Ala) 15, and (Ala) 20- As 
the number N of residues increases, the peak temperature and peak height also increase, 
and the peaks become more pronounced, suggesting the divergence of the specific heat 
in the thermodynamic limit (N —> 00). This behavior would be expected for a phase 
transition. We do not expect similar behavior for the other homo-oligomers, but did not 
study this question because of limited computer time. 

Finally, we investigate how each of the energy terms in Eqs. (16) - (20) varies as 
T changes. In Figure 6 we show each energy term as a function of temperature for 
(Ala) 10, (Val)io, (Gly)io, and (Ala) 20- In all cases, each term monotonically increases as 
T increases. The changes, however, are very small except for the Lennard- Jones term, 
Elj, indicating that E^j is the key factor for the folding of these homo-oligomers. ||1 1|| For 
(Ala) at (A = 10 and 20) we observe a steep change of E tot and Elj (and other terms) 
at the transition temperatures (T = 400 ~ 500 K), and the slope becomes steeper as N 
increases. These facts are reflected in the pronounced peaks in the specific heat in Figure 
5 and again suggest the existence of a phase transition for this polymer. Work is under 
progress to study the nature of this transition in more detail. 
Thermodynamics of Helix-Coil Transitions 

In Figure 7 we show the average helicity as a function of temperature for (Ala) 10, 
(Val)io, and (Gly)i . These were again calculated by the reweighting techniques of Eqs. (4) 
and (5). Here, < n > stands for < ur >. From the regular multicanonical production 
run, it was found that right-handed and left-handed helicities are essentially equal for 
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(Gly)io: < tir til >. The average helicity tends to decrease monotonically as the 

temperature increases because of the increased thermal fluctuations. Around the room 
temperature, (Ala)io is substantially helical (~ 75 % helicity), (Val)io is slightly helical 
(~ 30 % helicity), and (Gly)io is hardly helical (< 10 % helicity). This is consistent with 
the fact that Ala is a helix former and Gly is a helix breaker, while Val comes in between 
the two. For reasons described above we could not reweight our < n > data to T = 
K and we did not dare to extrapolate them from the safe temperature range to T = 
K. However, the qualitative behavior of < n > for decreasing temperature supports our 
conjecture that the ground state of both (Ala)io and (Val)io is an ideal helix, while this 
is not the case for (Gly) 10 . 

Again we investigated for (Ala) at the dependence of this quantity on the number of 
residues N. In Figure 8 we show as a function of temperature for (Ala)io, (Ala) 15, and 
(Ala) 20- The data indicate that the longer the homo-oligomer, the more helix-forming it 
is, as we guessed above. The difference between the residue numbers 15 and 20 is smaller 
than that with 10, indicating that the N — > 00 limit for (Ala) at may be reached already 
near N = 20. 

According to the Zimm-Bragg model, the average number of helical residues < n > 

and the average length < £ > of a helical segment are given for large N by 

< n > 1 1 - s 

N 2 2y /(l -sf + Asa ' 

< t > = 1 + , 23 , (23) 

1 - s + ^J(l - s) 2 + Asa 

where N is the number of residues, and s and a are the helix propagation parameter and 
nucleation parameter, respectively. Note that s > 1 implies that > \ (more than 50 
% helicity). From these equations with the values of and < t > calculated from the 
multicanonical production runs, one can obtain estimates of s and a parameters. The 
values at six temperature values T = 250, 300, 350, 400, 450, and 500 K for (Ala)io, 
(Ala)i 5 , (Ala) 2 o, (Val) 10 , and (Gly)io are listed in Table 4. The s parameters monoton- 
ically decrease as the temperature increases because the increased thermal fluctuations 
will decrease helical content as the temperature increases (see Figure 7). Note that the 
difference between N = 15 and 20 for (Ala) at is again small, suggesting that N = 20 is al- 
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(22) 



ready large enough. The difference between N = 10 and 20 is also within the errors. The 
s parameters at experimentally relevant temperature (~ 270 K) are about 1.6, 0.6, 0.2 for 
Ala, Val, Gly, respectively. These are in good agreement with the experimental results, || 
where they give s(Ala) = 1.5 ~ 2.19, s(Val)= 0.20 ~ 0.93 and s(Gly) = 0.02 ~ 0.57. The 
u parameters in Table 4, on the other hand, tend to be order of magnitude larger than the 
commonly assumed values (~ 10~ 3 ). Since the data for Ala have smaller errors and are 
more reliable than those for Val and Gly, we show the s and a values for Ala as functions 
of temperature in Figure 9. As discussed above, (Ala)io has a clear phase-transition sig- 
nal with a pronounced peak in specific heat at the transition temperature T c of ~ 430 K 
(see Figure 4). Above the transition temperature T c , (Ala) 10 is in the random-coil phase, 
and below that temperature it is in the helix phase (see also Figure 7). This transition 
temperature T c can also be identified with as the temperature where s — 1 holds (i.e., 
50 % helicity) in Figure 9a, which is about 410 K. The small disagreement is due to the 
arbitrariness in our definition of a helical state. As is clear from Figure 9b, in the helix 
phase (T < T c ) the a parameter for Ala is small and constant, but in the random-coil 
phase (T > T c ) o starts to grow as temperature increases. This growth of a values reflects 
the increased thermal fluctuations that prevent the formation of a long helix. That is, 
below T c cooperativity for helix-formation wins over thermal fluctuations but above T c 
thermal fluctuations win and no long helices can be formed. 

Figure 9 implies that as the number of residues N increases for (Ala) at, s values at 
a fixed temperature tend to increase and a values tend to decrease. Eqs. (22) and (23) 
are exact in the limit N — > oo. Hence, the s and a should be extrapolated in this limit. 
However, we found that the differences in the value of the s parameter at experimentally 
relevant temperature, around 273 K, are within the errors. Therefore we conjecture that 
iV = 10 is already large enough for determining this quantity. We expect that systematic 
errors due to our poor approximation of solvent effects are larger than the finite-size 
effects. The situation may be different for a. This quantity is decreasing as N increases 
and is < 0.06 for iV = 20, which is a more acceptable value. Obviously the reliability of 
our values for this quantity is much more limited by the small number of residues in our 
simulations. 
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Another way to investigate the limitations of our results due to finite number of 
residues is to look into the end effects of the homo-oligomers. Our analyses so far have 
been neglecting them. As for helical content, one expects to see fraying at the edges. 
In Figure 10 we show the percent helicity as a function of the residue number at 
T = 250 and 350 K for (Ala)i , (Val)i , and (Gly)io- Again these results were calculated 
by the reweighting techniques of Eqs. (4) and (5). We do observe fraying for all cases. 
The contrast is most outstanding for (Ala)io because it has high helicity. The increase of 
fraying as the temperature is raised is clearly seen for (Ala)i . 

The helix-coil transition can be further studied by calculating the free energy differ- 
ences AG = G H — G c , enthalpy differences AH, and entropy differences TAS between 
helix (H) and coil (C) states. Here, a conformation is considered to be in the helix state 
if it has a segment with helix length £ > 3. The free energy differences were calculated 
from 

AG=-RT\n <J ^ H> , (24) 
< N c > v 1 

where < Nh > and < Nc > are the average numbers of conformations in helix state and 
in coil state, respectively. The enthalpy differences were estimated from 

AH =< E H > - < E c >, (25) 

where < E H > and < E c > are average total potential energies < E tot > (see Eq. (16)) 
in helix and coil states, respectively. Finally, the entropy differences were obtained from 
AG and AH by the relation 

TAS* = AH - AG . (26) 

In Table 5 we list AG, AH, and TAS* at six temperatures T = 250, 300, 350, 400, 
450, and 500 K for (Ala)io, (Ala)i5, (Ala)2o, (Val)io, and (Gly)io- AG is negative for 
(Ala) io if T < 450 K. Hence, (Ala)i favors a helix state for T < 450 K. Again this 
limiting temperature is around the transition temperature T c of the phase transition 
from the helix phase to the random-coil phase discussed above. For (Val)i , a helix 
state is slightly favored if T < 300 K, and for (Gly)io, a coil state is favored at all 
temperatures. The reason for the variance in helix-propensities becomes clear by looking 
at the enthalpy differences AH for the different oligomers. For (Ala)i we observe large 
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enthalpy differences between helix and coil states which favor the helix state and win over 
the entropic term of opposite sign. On the other hand, the helix state is only slightly 
energetically favored for (Val)io and not for (Gly)io- For (Val)io enthalpy and entropy 
differences are of the same order, while for (Gly)io the entropic term wins and coil states 
are favored. We remark that the absolute values of TAS become larger as T increases 
because of the increased thermal fluctuations. 
CONCLUSIONS 

In the present work we have demonstrated the effectiveness of multicanonical algorithms 
and presented various quantities one can calculate by this method. This was done by 
taking the example of helix-coil transitions of amino-acid homo-oligomers. Our results 
show that we are able to calculate the helix propensities of amino acids and compare 
them with experiments. However, our results have to be taken as preliminary. They 
are both hampered by the small number of residues and even more the neglecting of 
solvent effects. We expect to obtain an even better agreement with experiments, once we 
incorporate more realistic energy functions with solvent effects included. 
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Table Captions: 

Table 1: Energies E (in kcal/mol) of the lowest-energy conformations in helix state and 
coil state obtained during each run of multicanonical annealing simulation. R-Helix and 
L-Helix stand for right-handed a-helix state and left-handed a-helix state, respectively. 
Helix length £ is also given for the conformations in helix state. The null entry implies 
that there was no conformation in that state with E < 35 kcal/mol. MIHC stands for 
minimized idealized helical conformation. 

Table 2: Dihedral angles (in degrees) of the lowest-energy conformations in helix state 
and coil state obtained in 10 runs of multicanonical annealing simulation. R-Helix and L- 
Helix stand for right-handed a-helix state and left-handed a-helix state, respectively. The 
symbols * and # indicate that the corresponding residue is respectively in right-handed 
and left-handed helix configuration. 

Table 3: Dihedral angles (in degrees) of the lowest-energy conformations in helix state 
and coil state obtained in multicanonical production runs. R-Helix and L-Helix stand for 
right-handed a-helix state and left-handed a-helix state, respectively. The symbols * and 
# indicate that the corresponding residue is respectively in right-handed and left-handed 
helix configuration. 

Table 4: Average number of helical residues < n >, average length of a helical segment 
< £ >, and the Zimm-Bragg s and a parameters as functions of temperature T (K). The 
numbers in parentheses represent errors. 

Table 5: Free energy differences AG(= Gh — Gc), enthalpy differences AH, and entropy 
differences TAS (all in kcal/mol) between helix (H) and coil (C) states as functions of 
temperature T (K). The numbers in parentheses represent errors. 
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Table 1. 





(Ala) io 


(Val)io 




(Gly)io 


State 


R-Helix 


Coil 


R-Helix 


Coil 


R-Helix 


L-Helix 


Coil 




E 




E 


E 




E 


E 


^ IX 


E 




E 


Run 
























1 


-9.1 


8 


3.7 


31.6 


3 


19.8 


4.3 


3 






-5.3 


2 


15.9 


4 


2.5 


6.6 


3 


5.9 


7.4 


3 


5.1 


3 


-4.8 


3 


-8.9 


8 


6.0 


12.9 


3 


14.8 


6.2 


3 


7.1 


5 


-5.3 


4 


-9.1 


8 


11.8 


5.3 


3 


6.1 


3.2 


3 


16.2 


3 


-5.5 


5 


-8.3 


8 


9.8 






5.9 


8.2 


3 


24.2 


3 


-4.6 


6 






-1.5 


4.6 


4 


6.4 


0.1 


3 


12.5 


3 


-7.4 


7 


0.2 


3 


-0.9 


14.4 


3 


8.7 


7.0 


3 


2.3 


3 


-3.6 


8 


-1.8 


3 


1.5 


4.3 


3 


5.1 


5.5 


3 


-1.3 


5 


-8.0 


9 


-9.1 


8 


3.8 


2.6 


6 


7.3 


9.9 


3 


5.5 


3 


-6.7 


10 


-8.2 


8 


7.0 


7.6 


3 


4.0 


27.3 


3 


-6.0 


6 


-5.4 


MIHC 


-8.8 8 




-3.4 8 




-3.9 8 


-3.8 8 





23 



(Ala) io 



Run 


State 


E (kcal/mol) 


9 


R-Helix 




9.1 


Residue 






X 




1 






-61 




2* 


-68 


-32 


-58 




3* 


-69 


-39 


-62 




4* 


-71 


-35 


-175 




5* 


-71 


-36 


180 




6* 


-73 


-34 


-179 




7* 


-71 


-36 


-54 




8* 


-72 


-37 


-62 




9* 


-71 


-40 


61 




10 


-145 


100 


58 








riR 


i L 


riL 




8 


8 









(Ala) io 



Run 


State 


E (kcal/mol) 


6 


Coil 




1.5 


Residue 


4> 


y, 


X 




1 


120 


-58 


57 




2 


-83 


72 


-178 




3 


-65 


110 


-64 




4 


-95 


-34 


-58 




5 


-153 


176 


63 




6* 


-57 


-49 


-59 




7* 


-71 


-50 


-60 




8 


-83 


74 


-178 




9* 


-82 


-22 


65 




10* 


-68 


-42 


177 






£r 


riR 


i L 


til 




2, 2 


4 









24 



Table 2. (continued) 
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Table 2. (continued) 
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Table 2. (continued) 
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Table 3. (continued) 
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Table 3. (continued) 
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R-Helix -5.2 



Residue 




y, 


1 


115 


-74 


2* 


-68 


-36 


3* 


-68 


-32 


4* 


-71 


-45 


5* 


-65 


-32 


6* 


-78 


-56 


7 


159 


68 


8 


-174 


-106 


9 


-164 


-86 


10 


-70 


-81 



£r n R i L n L 
5 5 



(Gly)io 





State 




E (kcal/mol) 




L-Helix 




-6.7 


Residue (j) 


V' 




1 


-22 


-53 




2 


-61 


95 




3# 


65 


34 




4# 


76 


35 




5# 


70 


45 




6# 


75 


39 




7# 


69 


31 




8* 


-87 


-38 




9 


147 


-33 




10 


-161 


-50 






Ir 


nR 


P-L n L 




1 


1 


5 5 
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Table 3. (continued) 



(Gly)io 

State E (kcal/mol) 

Coil -7.6 



Residue 




y, 


1 


-9 


66 


2 


57 


-106 


3* 


-72 


-28 


4* 


-71 


-30 


5 


-85 


-60 


6 


-179 


68 


7 


75 


-68 


8 


-130 


37 


9 


-81 


78 


10 


79 


74 



wr @L n L 
2 2 
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Table 4. 



Peptide 


T 


< n > 


< I > 




s 


a 


(Ala) 10 


250 


7.8(0.4) 


7.7(0.7) 


1.6 


[0.1) 


0.14(0.01) 




300 


7.5(0.6) 


7.2(1.0) 


1.5 


;o.2) 


0.13(0.02) 




350 


6.8(0.8) 


6.3(1.2) 


1.3 


;o.2) 


0.13(0.02) 




400 


5.5(0.9) 


4.7(1.0) 


1.1 


;o.2) 


0.14(0.02) 




450 


3.5(0.7) 


2.8(0.7) 


0.76 


;o.i2) 


0.18(0.05) 




500 


2.1(0.5) 


1.6(0.4) 


0.49 


;o.io) 


0.27(0.07) 


(Ala) 15 


250 


12.8(0.1) 


12.7(0.2) 


1.7 


;o.i) 


0.077(0.01) 




300 


12.6(0.2) 


12.2(0.4) 


1.6 


;o.i) 


0.072(0.01) 




350 


12.2(0.3) 


11.3(0.8) 


1.5 


;o.i) 


0.072(0.02) 




400 


11.4(0.4) 


9.8(1.0) 


1.4 


;o.i) 


0.076(0.06) 




450 


8.6(1.2) 


6.4(1.1) 


1.1 


;o.i) 


0.092(0.01) 




500 


5.5(2.1) 


3.5(1.7) 


0.77 


[0.24) 


0.17(0.08) 


(Ala) 20 


250 


17.5(0.3) 


15.9(0.9) 


1.7 


;o.i) 


0.051(0.01) 




300 


17.2(0.4) 


14.9(1.4) 


1.6 


;o.i) 


0.051(0.01) 




350 


16.7(0.5) 


13.2(1.5) 


1.5 


;o.i) 


0.052(0.01) 




400 


15.9(0.5) 


11.1(1.5) 


1.4 


;o.i) 


0.053(0.01) 




450 


14.6(0.7) 


8.8(1.2) 


1.3 


;o.i) 


0.058(0.01) 




500 


9.7(1.0) 


4.3(0.2) 


0.98 


(0.1) 


0.091(0.02) 


(Val)io 


250 


3.9(1.3) 


2.3(0.5) 


0.73 


[0-25) 


0.47(0.12) 




300 


3.4(1.2) 


2.0(0.5) 


0.60 


[0.22) 


0.57(0.30) 




350 


2.8(1.0) 


1.7(0.4) 


0.45 


[0.17) 


0.72(0.31) 




400 


2.3(0.9) 


1.4(0.5) 


0.36 


[0.18) 


0.85(0.40) 




450 


2.0(0.8) 


1.3(0.4) 


0.30 


[0.17) 


0.89(0.39) 




500 


1.7(0.8) 


1.1(0.4) 


0.26 


[0.16) 


0.89(0.36) 


(Gly)io 


250 


0.84(0.10) 


1.3(0.1) 


0.24 


[0.03) 


0.31(0.7) 




300 


0.81(0.03) 


1.2(0.1) 


0.20 


[0.14) 


0.38(0.4) 




350 


0.76(0.07) 


1.2(0.1) 


0.16 


[0.09) 


0.46(0.27) 




400 


0.68(0.09) 


1.1(0.1) 


0.14 


[0.05) 


0.48(0.14) 




450 


0.59(0.07) 


1.1(0.1) 


0.12 


[0.02) 


0.47(0.09) 




500 


0.52(0.05) 


1.1(0.1) 


0.11 


[0.02) 


0.45(0.05) 
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Peptide 


T 


AG 


AH 


TAS 


(Ala) 10 


250 


-3.9(0.8) 


-8.1(2.1) 


-4.2(2.1) 




300 


-3.1(0.7) 


-7.6(2.0) 


-4.5(2.0) 




350 


-2.3(0.6) 


-9.0(1.4) 


-6.7(1.4) 




400 


-1.2(0.6) 


-11.2(1.2) 


-10.0(1.4) 




450 


0.1(0.6) 


-11.9(1.5) 


-12.0(1.6) 




500 


1.4(0.6) 


-11.2(1.3) 


-12.6(1.7) 


(Ala) 15 


250 


-11.9(1.0) 


-25.2(2.2) 


-13.3(2.4) 




300 


-9.4(0.8) 


-23.6(2.7) 


-14.2(2.1) 




350 


-7.0(0.7) 


-24.1(2.4) 


-17.1(2.4) 




400 


-4.5(0.9) 


-24.3(3.2) 


-19.8(3.3) 




450 


-2.2(1.3) 


-19.9(4.5) 


-17.7(4.7) 




500 


-0.7(1.5) 


-13.3(6.4) 


-12.6(6.6) 


(Ala) 20 


250 


-22.1(0.6) 


-47.9(2.3) 


-25.8(2.3) 




300 


-17.2(0.6) 


-44.8(2.7) 


-27.6(2.7) 




350 


-12.9(0.4) 


-42.1(2.9) 


-29.2(2.9) 




400 


-8.8(0.3) 


-40.4(2.2) 


-31.6(2.2) 




450 


-4.9(0.4) 


-39.4(0.8) 


-34.5(0.9) 




500 


-1.5(0.3) 


-29.3(3.3) 


-27.8(3.3) 


(Val)io 


250 


-0.4(0.2) 


-2.1(0.7) 


-1.7(0.3) 




300 


0.0(0.5) 


-3.0(0.7) 


-3.0(0.7) 




350 


0.6(0.9) 


-4.1(1.5) 


-4.7(1.6) 




400 


1.3(1.1) 


-4.4(2.1) 


-5.7(2.3) 




450 


2.1(1.5) 


-4.4(2.5) 


-6.5(3.8) 




500 


2.8(2.0) 


-4.6(3.4) 


-7.4(5.4) 


(Gly)io 


250 


2.0(1.0) 


0.3(1.5) 


-1.7(1.8) 




300 


2.4(0.8) 


-0.3(1.1) 


-2.7(1.4) 




350 


2.9(0.8) 


-1.1(0.9) 


-4.0(1.1) 




400 


3.6(0.8) 


-1.8(1.6) 


-5.4(1.9) 




450 


4.3(0.9) 


-2.4(2.0) 


-6.7(2.6) 




500 


5.0(1.0) 


-3.0(1.8) 


-8.0(2.6) 
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Figure Captions: 

Fig. 1: Conformations with the lowest potential energy for each homo-oligomer obtained 
from multicanonical annealing runs and one regular multicanonical production run with 
200,000 MC sweeps. Stereoscopic views of the backbone structure are shown for (Ala)i 
(a), (Val)io (b), (Gly)i (c), and for the lowest-energy left-handed a-helix obtained for 
(Gly)io (d). In (a) and (b) residues from 2 to 8 and from 2 to 6 are respectively helical. 
The conformation in (c) is in the coil state. Residues from 3 to 7 are helical in (d). 
Fig. 2: Probability distribution P(E) of the energy as a function of temperature T for 
(Ala) io (a), (Val)io (b), and (Gly)i (c) obtained by one regular multicanonical production 
run with 200,000 MC sweeps. 

Fig. 3: Probability distribution P(E) of the energy for T =270 K as obtained by a canon- 
ical simulation (can) of 200,000 MC sweeps and a multicanonical simulation (mul) of 
200,000 MC sweeps for (Ala)i (a) and (Ala) 20 (b). 

Fig. 4: Specific heat as a function of temperature for (Ala)i , (Val)i , and (Gly)i . The 
values were calculated from one regular multicanonical production run of 200,000 MC 
sweeps. 

Fig. 5: Specific heat as a function of temperature for (Ala) at with iV = 10, 15, and 20. 
The values were calculated from one regular multicanonical production run of 200,000 
MC sweeps. 

Fig. 6: Average total energy E tot = E c + E HB + E LJ + E tor (C>) and averages of its com- 
ponent terms, Coulomb energy Ec (+), hydrogen-bond energy Ehb ( n ), Lennard- Jones 
energy E LJ (X), and torsion energy E tor (A) as a function of temperature T. The results 
are separately shown for (Ala)io (a), (Val)io (b), (Gly)io (c), and (Ala)2o (d). 
Fig. 7: Average helicity as a function of temperature T for (Ala)i , (Val)i , and 
(Gly)io- All values were calculated from one regular multicanonical production run of 
200,000 MC. 

Fig. 8: Average helicity as a function of temperature T for (Ala)^v with iV = 10, 15, 
and 20. The values were calculated from one regular multicanonical production run of 
200,000 MC sweeps. 

Fig. 9: Helix propagation parameter s (a) and nucleation parameter a (b) of the Zimm- 
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Bragg model as a function of temperature T for (Ala) ^ with iV = 10, 15, and 20. The 
values were calculated from one regular multicanonical production run of 200,000 MC 
sweeps. 

Fig. 10: Percent helicity as a function of residue number for (Ala)i , (Val)io, and (Gly)io 
at T = 250 K (a) and at T = 350 K (b). The values were calculated from one regular 
multicanonical production run of 200,000 MC sweeps. 
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